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^ Abstract. We examine the effect of subdividing the potential barrier along the reaction 

O coordinate on Kramers' escape rate for a model potential. Using the known supersymmet- 

^ I ric potential approach, we show the existence of an optimal number of subdivisions that 

\Q ' maximises the rate. 
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O . 1. Introduction 
O ■ 

. !-H ! The problem of surmounting a potential or, more generally, a free energy barrier is a classical 
^ ■ problem that appears in all processes having thermally activated kinetics. This problem was 
originally addressed by Kramers in the context of a bistable potential energy curve [1]. He 
provided an approximate solution for the rate of escape over the barrier in a high barrier 
low noise hmit. In the commonly encountered high friction limit, the bistable potential is 
usually parameterized in terms of the height of the barrier at the potential maximum and 
the width of, or the distance under, the barrier connecting the initial and the final states 
(potential minima). Since Kramers' original work, there has been a number of refinements 
as well as varied novel applications of his solution, and a large volume of literature exists on 
this [2,3]. 

There are, however, situations where the initial and the final states are separated by a 
barrier which is so high that the estimated reaction rate is very small, and yet the reaction 
actually turns out to proceed at a substantially higher rate. The enhancement is attributed 
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to the catalytic action, notably of an enzyme in a biochemical reaction that forms a 'tran- 
sition state' complex with the substrate giving a reduced barrier height [4]. We, however, 
envisage here an alternative scenario where the enzyme effectively reduces the activation 
energy by subdividing the reaction path into a number of discrete steps each requiring a 
much smaller barrier crossing. These subdivisions are expected to correspond to the discrete 
conformational/configurational changes of the macromolecules, protiens say. Besides look- 
ing for a physical consequence of the barrier subdivision, in its own right, the problem can 
be viewed as an excercise on rate processes in dissipative systems. In the present work we 
have considered the effect of the barrier subdivision on the reaction rate in the high friction 
limit. This we have done for a W-shaped model potential barrier whose subdivision can be 
well parameterized. Our analysis of the problem is based on the supersymmetric potential 
technique [5-7]. 

2. The Methodology 

It is sufficient for our purpose to note that in the high friction limit, Kramers' escape problem 
is one of solving the Smoluchowski equation (SE): 

where P(x,t) is the probability density associated with the particle position, U'{x) = ^ 
with U{x) being the 'double well'-potential, D is the diffusion constant and /3 = {kT)~^ is 
the inverse temperature. With the ansatz 

P(x,i) = 0(x)e-^^(^)/2g-At (2) 

the SE is converted to a Euclidean Schroedinger equation for 0: 

H+^+^E+4>+ (3) 
with if+ = A^A being positive semi-definite, where — X/D and. 

This Hamiltonian corresponds to the motion of a particle in the potential 

V^{x)^{^IpU'{x))' -IpU"{x). (6) 

For a high barrier, the escape rate is determined by the smallest nonzero eigenvalue, Ai = 
DE^, of the SE where E^ is the eigenvalue of the first excited state of Eq.(3). On the other 
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hand, this eigenstate is degenerate with the ground state 0^ of the 'supersymmetric partner 
potential' V^{x) given by 

v4x)^(^Ipu'{x))\Ipu"{x) (7) 

so that 0° = 0° with if_ = AA'^ and = E\. The problem thus boils down to 
finding the ground state eigenvalue of this 'partner' potential. 



3. The model and its solution 



3.1 The model potential and parameterization of subdivision 

For the sake of simphcity we consider a symmetric W-potential. For a full characterisation 
of this potential we require two parameters, namely, the height Uq and the width under the 
potential 2Lq (see Fig. la). Wc now subdivide the barrier between the initial and final 
states into a scries of smaller connecting barriers of many steps (sec Fig. lb). In order to 
examine the effect of barrier subdivision on the reaction rate systematically, it is necessary to 
parameterize the subdivision in a physically, reasonable manner. Consider the step located 
between Xn and x^+i (Fig. Id). We choose C/i, U2 and the associated widths a, h (shown 
in the figure) such that ^ = This choice simplifies the calculation further. Note that 
Xn+i — Xn = a + b. If we have a total of N such equally spaced steps from the top of the 
barrier on either side, then Lq = Na + {N — 1)6, while Uq = NUi — {N — l)U2. We introduce 
a parameter p defined as 

(N-m {N - i)b 



The aim is, given Uq and Lq, to find the escape rate for different values of barrier subdivi- 
sion consistent with the high barrier limit, i.e., for various values of N. Such parameterization 
is physically reasonable as it not only keeps the barrier height and its width fixed, it also 
keeps the area under the barrier approximately constant as N is varied. Introducing a 
dimensionless potential u{x) = ^pU{x), the 'supersymmetric partner potential' to V+{x), 
nsimely,V-{x), is then given by 

V4x)^{u'{x)Y + u"{x) (9) 
which for the considered potential takes the form 

V_{x)=(^) - — E [S{x-xn)- S{x-xn-a)] (10) 

where Ui = \l3Ui. Note that the potential V_{x) is a series of attractive and repulsive delta- 
potentials superimposed over a constant potential (see Fig. Ic). Changing the variable x to 
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y — ^ leads to a new Hamiltonian h- given by 

f)2 +N 

h_=alH_^-— + {ulf-2ul [Ky-yn)-5{y-yn-a,)]. (11) 

'^y n=-N 

where ao = Na, u\ — Nui, Oi = 61 = jf-^ and Un = n{ai + bi). With this, h^(f)^_{y) — 
e_ 4>^{y), where e_ is a dimensionless quantity equal to a^E_. 

3.2 Solution 



We use transfer matrix method to find the ground state energy. The ground state wave 
function 0'^ is of the form Ae'^"^ + Be^^ peaked around the positions of the delta potentials. 
Consider one period of the potential, say, the interval between yn and yn+i- Assume the 
wave function of the form 

01 (y) = e-'^^^-^") + S„ e^iy-y^)^ (12) 

for the interval + oi < y < yn and the wave function of the form 

(t>2{y) = Cne-^^y-y--"^^ + L>„e*^(^-^"-"i) (13) 

for the interval yn < y < y« + with k = [{u\y — e_]5. By matching the wave function 
at y^, i.e., (piiy-n) = <p2{yn), and by integrating Eq.(ll) around y„ noting that there is a 
negative delta-potential of strength 2u'^, i.e., 0'i(yn) — 02(yn) — 2ii?0i(yn) = 0, we get 

[dJ - V -a e^^- (1 - a)e'=- ) [bJ ' [bJ 

relating the two pairs of amplitudes (« = -^)- Next, matching the wave function having 
amplitudes C„ and Dn with the wave function having amplitudes An+i and -Bn+i (found 
in the interval yn + 0,1 < y < y„+i) and integrating Eq.(ll) around y„ + Oi noting that a 
positive delta-potential of strength 2u'( is located there, we get 

The transfer matrix, T, relating the amplitudes An+i, Bn+i to the amplitudes An and Bn'- 



(t:)^na '''' 

will then be a product of the matrices Ti and T2, i.e., T = T2T1. The amplitudes just 
before the end of the N-th step An, Bn are related to the amplitudes ^40, Bq at the top left 
side of the barrier by a product of N of these transfer matrices: 
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Symmetry of the potential about y = imphes that the ground state wave function is 
symmetric, i.e., (j)^{—y) = (f)^{y). Using matching and integration at and around the origin 
where there is a negative delta-potential with the symmetry property of the wave function 
relates Aq and Bq: 

Bo = \^Ao. (18) 

Since we are concerned with a bound state solution, B]\f = 0. Using this and Eq.(^) in 
Eq.([T7|) enables us to get 

(l-a)(T^)2i + (l + a)(T^)22 = 0, (19) 

the lowest positive solution of which gives us the value of e_ when uo, p and are specified. 
The expressions for the matrix elements (T'^)2i and (T^)22 are given in the Appendix. 

4. Results and discussion 

Now we consider the solution of Eg. ([T9|) . The result could be better appreciated if we 
compare it with the corresponding escape rate for the original W-potential with no barrier 
subdivision (Fig. la). Applying the same technique as above the equation corresponding to 
Eq.([T9|) to be solved is 

ao + (1 - ao) e^^° = (20) 

where cto = = [^o ~ ^°]^ with uq = ^PUq. In this case, e° = L^E^. The inverse of 

DE^ is the time required to go from one minimum to the other in Fig. (la) and we define 
the corresponding escape rate as DE^ for the original 'W potential. Then the ratio, /at, of 
the escape rate, DE^, over the potential with a certain barrier subdivision to that of escape 
rate, DE^, over the original W-potential is given by 

/^-§ = (1 + P)'^- (21) 

We call this ratio, f^, as the enhancement factor. It may be worthwhile pointing out here 
that we have used the first passage time from one minimum to the other in the original 
potential as a scale factor. This is because the subdivided potential is rugged on the 'down 
hill part' as well, which could give rise to a considerably different transit time compared to 
the situation if only 'sliding down' on a smooth line were allowed. 

There are only two parameters in our model, namely Mq, which is the total barrier height 
and p, which essentially represents the steepness of the local barriers. We chose Uq (so as 
to be in the high barrier limit) holding p fixed and explored the enhancement factor, /^r, 
for various values of barrier subdivision, A^. Fig. 2 shows plots of J'n versus A^ for three 
different choices of uq with fixed p(=0.8). For this case, the enhancement factor at the 
optimal barrier subdivision. Nop, increases as Uq is increased reaching a value as high as 35 
for Uq = 12.0, while Nop remains constant (here 9) suggesting that the steepness of the local 



5 



barrier determines Nop. In Fig. 3, we have varied p for a fixed value of mo(=9.0). In this 
case, both Nop and the enhancement factor increase as p is increased. ( Note that in both 
the figures N starts at values larger than 2 for higher values of p. This is due to the fact 
that for these values of p, the enhancement factor is less than unity for low values of N and 
correspond to the situation where U2 > Ui.) 

We have verified that these trends are general. Thus, there is an optimal value of barrier 
subdivision, Nop, at which the escape rate takes a maximum value. The existence of Nop may 
be readily understood by a reference to the potential V^{x). The binding energy of this 
localized ground state of the individual negative energy delta-functions and its lowering due 
to mutual overlap of the neighbouring bound states (banding effect) are oppositely affected 
by TV. 

It may be worthwhile to mention here that in addition to changing the terrain (say steep- 
ness) of the intermediate barrier (connecting the initial and final states) by the subdivision, 
the outer barriers were also made to change their steepness accordingly for the sake of sim- 
plicity (see Figs, la and lb). Because of this increased steepness, they become more confining 
than with their original slope and, thus, give rise to an overestimation of the enhancement 
factor. We have verified this by retaining the original slope of the outer barriers. However, 
the main features remain the same. On the other hand, due to the monotonic increase of 
the optimal enhancement factor with the barrier height, its value could be much larger than 
the ones considered here for barrier heights that exist in chemical and physical processes. 

We remark that this problem can be viewed, approximately, as that of finding the mean 
first passage time of a biased random walk [8] . However, this would imply assigning values 
to the forward and the backward transition rates for the individual sub-barriers taken in 
isolation as input, i.e. assuming that the potentials to the left and right sides of each sub- 
barrier are totally confining, and then using these input values to calculate the global escape 
rate for the coupled sub-barriers. We have found that while this gives an optimal barrier 
subdivision for the escape rate consistent with the present result, the enhancement factor 
is considerably over-estimated by this random walk approach. The present SUSY-based 
calculation goes beyond this uncontrolled approximation. 

It would be interesting to examine and optimize the effect of an athermal (possibly 
colored) noise (the 'blow torch' of Landauer [9,10]) on one of the steps of our subdivided 
potential curve. This is under investigation. 

In conclusion, we have shown that the Kramers' rate for the escape over a given potential 
barrier, in the high barrier high friction limit, can be substantially enhanced by subdivid- 
ing the barrier optimally. This might provide an alternative scenario for certain activated 
processes where the measured escape rate is substantially higher than that anticipated. 
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Appendix A 

To find the elements of the matrix we decompose the transfer matrix T as a product of 
three matrices; i.e. 

T = RAL (Al) 

so that A is a diagonal matrix whose diagonal elements are the eigenvalues of T. The 
matrices L and R are, respectively, made up of the left- and right-eigenvectors of T such 
that L R = R L = I. With this decomposition, 

T^=RA^L (A2) 

whose two elements of our interest, (T-^)2i and (T-^)22 are expressed as 

and 

(T )22 ^ . (A4) 

X± are the eigenvalues of T given by 

A± = ^(rn+T22 ±Q) (A5) 

with Tij as the matrix elements of T and Q = [(Tn - ^22)^ + 4:TuT2i]^. 
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Figure Captions 

Figure 1: (a) The model potential, (b) Plot of the subdivided model potential: U{x) 
versus x, when N — 3. Note the change in the slope of the left- and right-confining walls 
from that of (a), (c) Plot of V^{x) versus x (not to scale), (d) A figure showing a step of 
the subdivided potential found between the intervals x„ and Xn+i- 

Figure 2: Plots of /n versus N for three different values of uq (i.e. 6.0, 9.0, 12.0) with 
fixed p{— 0.8). Note that all the optimal values occur &t N — 9. 

Figure 3: Plots of /jv versus for four different values of p (i.e. 0.4, 0.6, 0.8, 0.9) with 
fixed Uo{= 9.0). 



8 



40 



30 - 



Uo=12.0 



X 



X 



X 



o o 



X 



o 



o 



X 



o 



20 - 



o 



o + + + + 



10 



° 0^0=9-0 



X 



o 



o 



o 



o 



o 



o 



o 



o 



Uo=6.0 



o 



o 



5 



10 15 20 25 



N 



40 



30 - 



20 



10 - 








5 10 15 20 25 30 



N 



